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We address here a few classical lattice-spin models, involving n—component unit vectors (n = 2, 3), asso¬ 
ciated with a _D—dimensional lattice D = 1, 2, and interacting via a pair potential restricted to nearest 
neighbours and being isotropic in spin space, i.e. defined by a function of the scalar product between the in¬ 
teracting spins. When the potential involves a continuous function of the scalar product, the Mermin-Wagner 
theorem and its generalizations exclude orientational order at all finite temperatures in the thermodynamic limit, 
and exclude phase transitions at finite temperatures when D = 1; on the other hand, we have considered here 
some comparatively simple functions of the scalar product which are bounded from below, diverge to -l-oo for 
certain mutual orientations, and are continuous almost everywhere with integrable singularities. Exact solutions 
are presented for D — 1, showing absence of phase transitions and absence of orientational order at all finite 
temperatures in the thermodynamic limit; for D = 2, and in the absence of more stringent mathematical re¬ 
sults, extensive simulations carried out on some of them point to the absence of orientational order at all finite 
temperatures, and suggest the existence of a Berezinskii-Kosterlitz-Thouless transition. 

PACS numbers: 05.50.-l-q, 64.60.-i, 75.10.Hk 


I. INTRODUCTION 

The study of lattice spin models, both classical (on which 
we shall be concentrating here) and quantum, is an important 
chapter of Statistical Mechanics, where a number of mathe¬ 
matical results have been obtained, entailing absence or exis¬ 
tence, and sometimes type, of phase transitions at finite tem¬ 
peratures, depending on lattice dimension, number of spin 
components, range and symmetry of the interaction. 

The Mermin-Wagner theorem was first proven nearly 50 
years ago in a quantum setting, for the isotropic spin Heisen¬ 
berg model with finite-range exchange interactions [1], and 
later extended by various Authors in a number of directions, 
e.g. to the classical setting, to other functions of the scalar 
product, or to longer-ranged interactions [2-5]; see also a sub¬ 
sequent Review in Ref. [6]. 

In the classical case, the Mermin-Wagner theorem and 
its generalizations [7-10] hold for lattice-spin models, con¬ 
sisting of n— component unit vectors (n > 2), associated 
with a ZJ—dimensional lattice Ajy (D = 1,2 and typically 
Ajj = Z^), and interacting via pair potentials which are 
isotropic in spin space, and usually translationally invariant 
(on the other hand, mathematical results have also been ob¬ 
tained which do not need any translational invariance [11- 
16]); the distance dependence is usually taken to be suitably 
short-ranged. Their orientational dependences are defined by 
some functions of the scalar product between interacting spin 
pairs: the earlier mathematical results were obtained for rather 
smooth functions (simple polynomials), and conditions were 
later gradually relaxed, i.e. to the milder request of continuity, 
and, in some cases, even to less regular functions [9, 10]. 

More explicitly, continuity is required in Refs. [9, 10], 
and some singularities are also allowed for in Ref. [9]; we 
are restricting our present discussion to finite-range (actually, 
nearest-neighbour) interactions, and notice that mathemati¬ 


cal results are known for long-range interactions as well (see, 
e.g., Refs. [8-10], and others quoted therein). 

To fix notation and ideas, let Wj = (wj, w?, • • • 
denotes the n-component unit vector (spin) associated with 
the j—th lattice site, with dimensionless coordinate vector 
Xj £ Z^; two-component spins are parameterized by usual 
polar angles ipj, and three-component spins are parameter¬ 
ized by usual spherical angles {6j, fij). Here and in the fol¬ 
lowing the interaction will be restricted to nearest neighbours 
and defined by 

T> = $(r) = ei^(r), T = Tyfe = Wy • Wfe, (1) 

where e > 0 denotes a positive quantity setting energy and 
temperature scales (i.e. T = ksTK/^, where 7k denotes the 
temperature in degrees Kelvin), and to be scaled away from 
the following formulae. For 2—component spins, it will prove 
notationally convenient to define 

T = cos A, A = Ajk = T’j - (2) 

When F{t) is a continuous function of its argument, the 
above theorems entail absence of orientational order in the 
thermodynamic limit at all finite temperatures [9]; when 
D = n = 2, and under additional conditions, a Berezinskii- 
Kosterlitz-Thouless (BKT), or, in more general terms, a BKT- 
like transition can be proven to exist [17-24]; the term “BKT- 
like” is used here to indicate a transition to a disordered low- 
temperature phase possessing slowly decaying correlations re¬ 
sulting in infinite susceptibility; in thermodynamic terms, the 
transition may be of infinite order (as in the more common, 
originally studied BKT case [17-22]; it was also later proven 
[23, 24] that it can turn first-order under certain conditions. 

Cases where F{t) possesses some singularity have been 
studied far less extensively (see also below). In fact one can 
envisage a multitude of singular interactions: models involv¬ 
ing a finite number of jump discontinuities, as in sign or step 
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models, are discussed in Appendix A; another family, also dis¬ 
cussed there, involves constrained models, where whole re¬ 
gions of configuration space are excluded. We have chosen 
to start our investigation, so to speak, somewhere in between 
these two cases, from functional forms containing slowly di¬ 
vergent terms which do not disturb thermodynamics, i.e. from 
functional forms being bounded from below, continuous al¬ 
most everywhere, slowly diverging to -|-oo for one (or a few) 
mutual orientations, and possessing integrable singularities. 
Thus the present paper addresses a few models whose func¬ 
tional forms are defined by 

^(t) = — ln(l-I-r), n = 2, (3a) 

kk(T) = — ln(l-b t), n = 3, (3b) 


X{t) = -ln(|r|), 

n = 3. 

(3c) 

In due course, comparisons will also be made with their ex¬ 
tensively studied counterparts defined by 

F(t) = -r, 

n = 2, 

(4a) 

F(r) = -r, 

n = 3, 

(4b) 

F(r) = -P 2 (t), 

n = 3, 

(4c) 


respectively, and simply referred to as “regular counterparts”. 

Some models bearing similarities to ours [Eq. (3b)] have 
been investigated previously in the literature [25-29]. More 
recent studies showed that such classical models are effec¬ 
tive models obtained via mappings from quantum-mechanical 
treatments [26-28]. The above singular models [Eqs. (3)], as 
well as some generalizations and linear combinations of them, 
can be solved exactly when D = 1, allowing one to obtain 
thermodynamic and structural quantities in closed form; these 
are worked out in Appendix A, where other singular models, 
such as step or sign model and constrained ones are addressed 
as well. The three models in Eqs. (3) are studied by extensive 
Monte Carlo (MC) simulation for D = 2 so as to explore the 
thermodynamic behavior of these models, on the one hand, 
and to unveil potential effects of the singularities in compari¬ 
son with their regular counterparts, on the other hand. 

The rest of the paper is organized as follows: in Sec. II we 
further discuss the singular models; our simulation methodol¬ 
ogy for D = 2 is discussed in Section III along with with brief 
details on the finite-size approach we employ for the analysis 
of the simulation data. In Sec. IV we present the simulation 
results and finite-size scaling analysis used to extract the criti¬ 
cal behavior for the models under consideration. We conclude 
the paper with Sec. V where we summarize our results. 

II. REMARKS ON THE POTENTIAL MODELS 

Both V and W attain their minimum at Tmin = 1, and 
slowly diverge to -boo as t —— 1; (r) attains its minima at 

Tmin = ±1, and slowly diverges to -boo as r —>■ 0; the above 
functions are bounded from below, continuous almost every¬ 
where, and possess integrable singularities; in these cases, an 


interaction diverging to -boo is still compatible with the ther¬ 
modynamics and, by its very functional form, it can be ex¬ 
pected to enforce some strengthening of short-range correla¬ 
tions. On the other hand, changing the sign in front of the 
“In” from ” to “-b” in (any of) Eqs. (3) would produce a 
rather dramatic effect, i.e. it would cause a divergence to —oo 
for some mutual orientations, and hence make the modified 
model not well defined at low temperatures [25]. 

Series expansions of Eqs. (3) can be written down, i.e. 

V{t) = ln(2) -b lim V„, 

(J—^OO 

Vq = 2^ — cos{lAjk), 0 < Ajk < TT- (5a) 
1=1 

Wt) = lim W„, 

q—¥oo 

^ f-lV 

= - (5b) 

X{r) = lim A'g, 

q—^oo 
^ 1 

= o< (5c) 

1=1 

each Xq is a polynomial in | t |, where the coefficient in front 
of |rp bears the sign (—1)*; in other words sign alternation is 
a common feature of the three above expansions. Any of the 
above truncated expansions [Eqs. (5)] is a continuous function 
of r which, by the Mermin-Wagner theorem and its general¬ 
izations [9, 10], produces orientational disorder at all finite 
temperatures; let us now consider a generalization of Vq, i.e. 

g 

•^9 = E c/cos(/A), (6) 

1=1 

where c; denote arbitrary real coefficients; the Mermin- 
Wagner theorem can be applied here as well; moreover, for a 
general ferromagnetic interaction (where all the coefficients c; 
are < 0), one can prove BKT behavior, based on its existence 
for Eq. (4a) [17] and on correlation inequalities, and also ob¬ 
tain a rigorous lower bound on the BKT transition temperature 
(see Ref. [21] and others quoted therein); unfortunately, the 
alternating signs in Vq prevent us from using this approach in 
general. Let us also mention in passing a simple specific case 
of Eq. (6), defined by 

C /2 = Cl cos A-b C 2 cos(2A), Cl < 0, (7) 

where C 2 can both be negative or sweep a suitable range of 
positive values; the model was studied by various Authors in 
the Literature (see Refs. [30, 31] and others quoted therein), 
also in the equivalent version [32, 33] (recall Appendix B) 

Qi = C 2 cos(2A) -b C 4 cos(4A), C 2 < 0; (8) 

simulation or spinwave evidence of BKT behavior was ob¬ 
tained in various cases, and estimates of the BKT transition 
temperature obtained for cases where the above mathemati¬ 
cal treatment applies [30, 32] were later shown to agree with 
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the named lower bound [21]. It proves convenient to compare 
each singular interaction potential [Eqs. (3)] with its regu¬ 
lar counterpart [Eqs. (4)], and with some truncated expansion 
[Eqs. (5)]; this is done in EIGs. 1, 2 and 3. These are found to 
exhibit a common feature; on the one hand, the singular inter¬ 
actions diverge rather slowly for appropriate mutual orienta¬ 
tions; on the other hand, in a broad minimum-energy region, 
the growth of the singular interaction energy as t moves away 
from the corresponding Tmin is recognizably slower than for 
its regular counterpart, and then it becomes faster and faster 
outside this region; the changeover takes place about r = 0 
(V and W model), or r Ri i {X model); a somewhat simi¬ 
lar behavior can also be seen for some (convergent) truncated 
expansions, and seems to reflect the above sign alternation. 



FIG. 1. (Color online) Comparison beteween the singular model V, 
its regular counterpart, and some truncated expansions [Eqs. (3a), 
(4a), (5a)], as functions of the angle A between the two spins. Mean¬ 
ing of symbols: red continuous line: regular counterpart; blue dashed 
line: model E; magenta dotted line: V 4 ; hrown dash-dotted line: 
Vio. 

What happens when the underlying lattice is taken to be 2- 
dimensional? The functional forms under investigation here 
[Eqs. (3)] diverge to -Foo for some mutual orientations, and, 
on the other hand. Refs. [9, 10] address the general case of 
continuous functions of the scalar product and Ref. [9] can 
even allow for some singularities; as far as we could check, 
the divergent behavior of the models under investigation here 
does not fit into the framework of weak singularity conditions 
used in section 2.2 of Ref. [9]. More explicitly, based on 
the series expansion in Eq. (5a), one could try to realize a 
decomposition of V (r) along the lines of Ref. [9], (sect. 2.2, 
around their Eqs. (24) to (26), page 441^43), by choosing a 
(large) positive integer q and rewriting Eq. (3a) as 

E(r) =ln(2)-f V, + r,; (9) 

the divergent term would then be positive around A = tt, 
and its sign would not agree with the hypotheses stipulated for 
theorem 1, singular case, in Ref. [9], where the small singular 
term in the interaction is written (their notation) 

-v{(j)), v{-) > 0 . 



FIG. 2. (Color online) Comparison beteween the singular model W, 
its regular counterpart, and some truncated expansions [Eqs. (3b), 
(4b), (5b)], as functions of the scalar product r between the two 
spins. Meaning of symbols: red continuous line: regular counter¬ 
part; blue dashed line: model IT; magenta dotted line: W 2 ; brown 
dassh-dotted line: VV 4 . 
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FIG. 3. (Color online) Comparison beteween the singular model X, 
its regular counterpart, and some truncated expansions [Eqs. (3c), 
(4c), (5c)], as functions of the scalar product r between the two spins. 
Meaning of symbols: red continuous line: regular counterpart; hlue 
dashed line: model X\ magenta dotted line: X 2 ; hrown dash-dotted 
line: X 4 . Notice that the quantity —1 has been added to X as well as 
to the two truncated expansions, in order to ease comparison. 


Thus there appears to be no available mathematical theorem 
entailing a Mermin-Wagner-type result in this case, although 
it has been conjectured (expectation is not calculation) that, in 
the thermodynamic limit, orientational order is also destroyed 
at all finite temperatures; (see. e.g. Ref. 13 in Ref. [28]); 
on the other hand, at least for the V case, one might expect 
a BKT behavior, since the singularity of the potential should 
ultimately strengthen short-range correlations. 
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III. SIMULATION ASPECTS AND FINITE-SIZE SCALING 
THEORY 


For D = 2, the three models V, W and X [Eqs. (3)] were 
treated by simulation. Calculations were carried out using 
periodic boundary conditions, and on samples consisting of 
N = particles, with L = 40,60, 80,100,120,160. Sim¬ 
ulations, based on standard Metropolis updating algorithm, 
were carried out in cascade, in order of increasing temperature 
T; equilibration runs took between 25000 and 50000 cycles, 
where one cycle corresponds to 2N attempted Monte Carlo 
steps, including sublattice sweeps (checkerboard decomposi¬ 
tion [34—37]), and production runs took between 500000 and 
1500000. 

Subaverages for evaluating statistical errors were calculated 
over macrosteps consisting of 1000 cycles. Calculated quan¬ 
tities include the potential energy (in units e per particle), and 
derivative with respect to temperature based on the fluctuation 
formula 


U* = 


W 

N '' 


( 10 ) 


and 


C* = ^ 


(11) 


with 


H= Y. Fir.k), ( 12 ) 

where J 2 {j<k} denotes sum over all distinct nearest- 
neighbouring pairs of lattice sites. 

As for orientational quantities, such as mean magnetization 
and corresponding susceptibilities [38, 39], they can be ex¬ 
pressed in general by 


Notice that Eq. (14a) involves a true ordering transition tem¬ 
perature Tc'. in our case, for models V and W, we found con¬ 
sistent evidence of the absence of orientational order at all 
finite temperatures (see also following Section), i.e. Tc = 0, 
and selected the definition of xi accordingly. Model X [Eq. 
(3c)], on the other hand, possesses even symmetry, and its 
second- and fourth-rank order parameters P 2 and P 4 , as well 
as the corresponding susceptibility X 2 , were calculated as dis¬ 
cussed in Ref. [40]; notice that, in this case 


X2 < PN. (15) 

We also calculated various short-range order parameters, de¬ 
fined by 


CTJ = {£j{Tjk)), (16) 

measuring correlations between corresponding pairs of unit 
vectors associated with nearest-neighbouring sites; here 
£j{t) denote appropriate orthogonal polynomials [see Eq. 
(A5) in Appedix A], and we chose J = 1, 2 for both V and 
W models, and J = 2,4 for the X model. 

In the quest for the possible occurrence of a phase tran¬ 
sition in the models investigated here, we will analyse the 
simulations data via the finite-size scaling (ESS) theory for 
continuous phase transitions - second order and BKT (infinite 
order) [37, 41, 42]. According to ESS hypothesis when a sys¬ 
tem is restricted to a finite geometry (a square of area in 
the present case) its thermodynamic quantities acquire a size 
dependence with a behavior that is tightly related to the order 
of the phase transition. It is worth mentioning that finite-size 
effects become important when the correlation length is of the 
same order as the linear size of the system. To be more specific 
we give details based on the behavior of the susceptibility. 

In the vicinity of a bulk critical point Tc the (magnetic) 
susceptibility diverges against the reduced temperature t = 
1 — ^ <C 1 according the scaling law xi ^ \t\~^ with the 
critical exponent 7 > 0. Eor a finite-size system it turns into 


N 

P = ^Wfc, (13a) 

k^l 

M=j^{\P\), (13b) 

M 2 = ^(P-P), (13c) 


r ^ (M 2 - iVM2) , T <Tc 
Xi = \ , (14a) 

[ /3M2, T>Tc 

where /3 = 1/T, and Tc denotes the critical temperature; since 
|P| < N [Eq. (13a)], we have 


Xi{L,T) = (17) 


where v measures the degree of divergence of the distance 
over which the spins are correlated, i.e. the correlation length 
f ^ with 1 / > 0. The function Oy^{x) is a universal 

function depending on the gross features of the system, but 
not of its microscopic details. 

On the other hand, when a BKT transition takes place, the 
susceptibility of the bulk system diverges exponentially 


Xbkt 


exp 


('x ( T ’ — Tbkt) 


1 

2 


Tbkt ^ T 


(18) 

as we approach Tbkt and is infinite in the BKT phase with 
a quasi-long range order. Eor a finite system however the di¬ 
vergence is rounded and the susceptibility is finite [Eqs. (14b) 
and (15)]. In the vicinity of the bulk BKT temperature the cor¬ 
relation length is proportional to the system’s linear size and 
the susceptibility scales like 


(14b) 


Xbkt 


M 2 < N and xi ^ P^- 


j^2 — r]BKT(T) 


(19) 
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At the transition temperature rjBKT = |- 

Expressions (17) and (19) are valid asymptotically in the 
vicinity of the transition temperature i.e. when both the sam¬ 
ple size L and the correlation length ^ are very large, but their 
ratio I; is finite. In this limit the universal scaling behavior is 
not affected by the finite-size effects. 


IV. SIMULATION RESULTS AND ESS ANALYSIS 

Simulation results obtained for the three investigated mod¬ 
els turned out to exhibit broad qualitative similarities, to be 
contrasted to their regular counterparts (see following discus¬ 
sion). 


A. The magnetic models V and W 

Simulation results for various observables, obtained for the 
two models V and W, were found to exhibit a recognizable 
qualitative similarity over a wide temperature range, so that, in 
some cases, only V results will be presented in the following. 

Simulation data for the potential energies of both models 
(not shown here) were found to evolve with tempereture in a 
gradual, monotonic way, and to be essentially independent of 
sample sizes, to within statistical errors falling below 0.1%. 



T 

FIG. 4. (Color online) The specific heat of model V for different 
sample sizes against temperature; statistical errors (not shown here) 
range between 1 and 5%. Meaning of symbols: red circles: L — 40; 
green squares: L = 60; blue triangles: L — 80; magenta diamonds: 
L = 100: cyan crosses, red asterisks: L = 160. 

As for the configurational specific heat C* (see FIG. 4 for 
model V, and FIG. 5, for model W), related to thermal fluctu¬ 
ations of the potential energy, the plots showed that C* starts 
with a maximum at T = 0, and first decreases to a broad 
minimum (say at T'); it then increases to another maximum 
(say at T"); here the associated statistical errors range be¬ 
tween 1 and 5%, and results are only mildly affected by sam¬ 
ple size. We found T' « 0.75, T" 1.2 for the V model. 


and T' k. 0.4, T" r; 0.62 for the W counterpart; upon ex¬ 
trapolating the low-temperature results to T = 0, we estimate 
the corresponding zero-temperature values to be ^ and 1, re¬ 
spectively; notice also that the zero-temperature value for the 
W model (but not for the V model) corresponds to the global 
maximum; on the other hand, T" for the V model (but not 
for the W model) corresponds to the global maximum. The 
same behaviour was found by estimating the specific heat via 
numerical differentiation of the internal energy. 



T 

FIG. 5. (Color online) The specific heat of model W for different 
sample sizes against temperature; statistical errors (not shown here) 
range between 1 and 5%; same meaning of symbols as in FIG. 4. 

A finite-size analysis of the configurational specific heat ac¬ 
cording to corresponding scaling behavior compatible with 
(17) ruled out the existence of a second order phase transi¬ 
tion in both models. A similar analysis was performed on the 
magnetization and the susceptibility for both models, but no 
scaling was achieved. 

Simulation results for the magnetization obtained with both 
models (see e.g. FIG. 6 for model V) showed a decreasing 
behavior as a function of temperature for a given sample size; 
at each examined nonzero temperature, they kept decreasing 
with increasing sample size; low-temperature results appear 
to extrapolate to M = 1 at T = 0 for all examined sample 
size, as expected. 

Low-temperature simulation results for M and for both 
models V and W were found to exhibit a power-law decay 
with increasing sample size; recall that the spin-wave analy¬ 
sis worked out in Ref. [43] for the regular counterpart [Eq. 
(4a)] predicts the low-temperature result 

_T_ 

M « {2L^) 8^ . (20) 

Our data at a given temperature were well fitted in a log-log 
scale by the relation 

InM = —alnL-F 6, a > 0, (21) 

where the ratio was found to increase with temperature, 
and to become constant in the low-temperature limit. 
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T 

FIG. 6. (Color online) Simulation results for the magnetisation of 
model V obtained with different sample sizes; same meaning of sym¬ 
bols as in FIG. 4. 


The thermal fluctuations of the magnetization for both mod¬ 
els V and W i.e. their magnetic susceptibilities (actually 
Inxi) are presented in FIGs. 7 and 8. At low tempera¬ 
tures the susceptibility keeps growing with sample size for 
both models, within the constraint of ( 14b), whereas at higher 
temperatures it becomes independent of sample size; the tem¬ 
peratures Tch where this change of scaling behavior first be¬ 
comes recognizable are Tch ~ 1.3 > T" for model V, and 
Tch ~ 0.56 < T" for model W, respectively. 

This specific behavior suggests a BKT transition from a 
quasi-long range ordered phase at low temperatures to a dis¬ 
ordered phase at higher ones. Assuming such a transitional 
behavior, we have fitted the data of the largest sample size 
(L = 160) to expression (18) for the bulk susceptibility and 
found the results of Table I, as crude estimates (see also be¬ 
low). 


TABLE I. Estimates of the parameters in Eq. (18) obtained by fitting 
to data for the largest sample size for models V, W and X assuming 
they exhibit a BKT transition. 


Model 

In ( ax ) 

K 

Tbkt 

V 

-3.21 ±0.15 

5.29 ±0.15 

0.873 ± 0.007 

W 

-3.93 ±0.14 

5.40 ±0.14 

0.259 ± 0.007 

X 

-2.33 ±0.05 

3.03 ± 0.04 

0.347 ± 0.003 


We analyzed the behavior of the susceptibility xi according 
to the finite-size scaling ansatz (19) in the vicinity of T = 0.9 
for model V and of T = 0.28 for model FF; we first car¬ 
ried out a linear fit of Inxi vs. liiT and estimated the crit¬ 
ical exponent p from the slope of the curves corresponding 
to different temperatures. The values obtained are presented 
in Tables II and III, for models V and W, respectively. A 
nonlinear fit, based on Eq. (19) was performed as well, and 
yielded results in agreement with these ones. Thus the transi¬ 
tion temperatures are most likely at Tbkt = 0.910 ± 0.005 



FIG. 7. (Color online) Simulation results for the susceptibility xi 
of model V obtained with different sample sizes; same meaning of 
synbols as in FIG. 4; Assuming a BKT transition and fitting the 
largest sample size L = 160 (upper continuous orange curve) to 
the bulk behavior of the susceptibility leads a transition at Tbkt = 
0.883 ± 0.007. 



T 


FIG. 8. (Color online) Simulation results for the susceptibility xi 
of model W obtained with different sample sizes; same meaning of 
symbols as in FIG. 4. 


and Tbkt = 0.275 ± 0.005 for models V and W, respec¬ 
tively. The discrepancy between these values and those in Ta¬ 
ble I points to the presence of huge finite-size effects: recall 
that Eq. (18) holds in the thermodynamic limit only, but was 
applied here to the largest investigated sample size in the hope 
to gain insights in the transitional behavior of the models con¬ 
sidered here. 

Eor the regular counterpart of model V the configurational 
specific heat was found to exhibit a sharp maximum at about 
15% [43] above the BKT transition. In Refs. [44,45] we have 
investigated the impact of diluted random impurities on the 
transition temperature. In Ref. [44] we have found a broad 
peaks about 5% above the BKT transition, and in Ref. [45] 
we found a sharper one about 2% above the transition temper- 
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ature. Here we find a maximum at about 40% above Tbkt- 
All these results show that the maximum of the specific heat is 
always above the transition temperature. As for Tch, we could 
not find in the Literature any estimate for the regular counter¬ 
part [Eq. (4a)]; thus additional simulations were run for the 
named regular model, carried out with the same sample sizes 
as for the three singular models, and using overrelaxation [46- 
50]; the estimate Tch ~ 1-05 was obtained. 


TABLE II. Estimates of rj for model V obtained via a log-log fit 
according to Eq. (19) for different temperatures along with the cor¬ 
responding error Srj. 


T 

0.890 

0.895 

0.900 

0.905 

0.910 

0.915 

0.920 

V 

0.246 

0.244 

0.248 

0.240 

0.250 

0.251 

0.257 

Sr) 

0.007 

0.006 

0.006 

0.006 

0.005 

0.004 

0.006 


TABLE III. Estimates of rj for model W obtained via a log-log fit 
according to Eq. (19) for different temperatures along with the cor¬ 
responding error Srj. 


T 

0.265 

0.270 

0.275 

0.280 

0.285 

0.290 

0.295 

V 

0.239 

0.243 

0.249 

0.257 

0.271 

0.273 

0.285 

Sr) 

0.004 

0.004 

0.004 

0.006 

0.006 

0.005 

0.007 



T 

FIG. 9. (Color online) Simulation results for the short-range order 
parameters (16) of model V obtained with the largest sample size; 
meaning of symbols: red diamonds (upper curve) tri; green circles 
(lower curve): < 72 . 

Simulation data for the short-range order parameters de¬ 
fined in ( 16) were found to be independent of sample size, and 
to decrease with temperature in a gradual and continuous way, 
paralleling the potential energy data; results obtained with the 
largest sample size of model V are collected on FIG. 9. 


B. The two-dimensional nematic model X 

Simulation results for the X model were also found to ex¬ 
hibit a remarkable qualitative similarity with the ones obtained 
for their magnetic counterparts. Data for the potential energy 
(not shown) as well as for the short-range order parameters 
(FIG. 10) were found to be independent of sample size, and to 
evolve with temperature in a gradual and monotonic way. The 
temperature dependence of the specific heat corresponded to 
its magnetic counterpart (FIG. 11); here also the associated 
statistical errors were found to range between 1 and 5%, and 
the results appeared to be only mildly affected by sample size. 
The plot started with the value 1 at T = 0, decreased with in¬ 
creasing temperature reaching a broad minimum at T' Ri 0.3, 
and then its global maximum at T" r; 0.5. it is worth men¬ 
tioning that a quite similar behavior was obtained by numer¬ 
ical differentiation of the potential energy. Notice also that, 
in the three cases, sample-size effects on the results become 
more pronounced about T". Here we anticipate that neither 
the results for the specific heat nor those corresponding to the 
second-rank order parameter P 2 or to the susceptibility X 2 
could obey the scaling behavior characteristic of a second or¬ 
der phase transition. 



0.0 0.2 0.4 0.6 0.8 

T 

FIG. 10. (Color online) Simulation results for the short-range order 
parameters (16) of model X obtained with the largest sample size; 
meaning of symbols: red diamonds (upper curve) G 2 \ green circles 
(lower curve): 0 - 4 . 

Simulation results for the order parameters Pj, (J = 2, 4) 
were also found to decrease with increasing temperature for 
each sample size, and to decrease with increasing sample size 
at each nonzero temperature (FIG. 12 and FIG. 13). At all 
investigated temperatures the results for the nematic order pa¬ 
rameters Pm, (M = 2,4) exhibited a power-law decay with 
increasing sample size. At a given temperature these were 
well fitted to the corresponding relations 

hiPj =-bjilnL + bjo, bji > 0. (22) 

The coefficients bji{T) were found to increase with T, and 
to become proportional to T to within statistical errors in the 
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T 

FIG. 11. (Color online) The specific heat of model X for different 
sample sizes against temperature; statistical errors (not shown here) 
range between 1 and 5%; same meaning of symbols as in FIG. 4. 


low temperature region. The results obtained from Eq. (22) 
show that both order parameters vanish in the thermodynamic 
limit i.e. L oo; such a behavior is in agreement with the 
spin wave theory for magnetic systems discussed above. 



0.0 0.2 0.4 0.6 0.8 

T 

FIG. 12. (Color online) Simulation results for the second-rank order 
parameter P 2 of model X obtained with different sample sizes; same 
meaning of symbols as in FIG. 4. 

Simulation results for In X 2 versus T (FIG. 14) showed a 
low-temperature regime where they kept increasing with in¬ 
creasing sample size, and then became independent of sam¬ 
ple size at higher temperatures; the temperature Tct where 
this change of scaling first becomes recognizable was T^h ~ 
0.45 < T"; this behavior also parallels the one observed for 
the two magnetic counterparts. 

By fitting the data obtained at high temperatures for our 
largest sample size (L = 160) to expression (18) of the sus¬ 
ceptibility, we obtain the results reported in Table I with a 
transition temperature Tbkt = 0.347 ± 0.003. 

Upon applying the finite-size-scaling analysis with data for 



0.0 0.2 0.4 0.6 0.8 

T 

FIG. 13. (Color online) Simulation results for the fourth-rank order 
parameter P 4 of model X obtained with different sample sizes; same 
meaning of symbols as in FIG. 4. 



FIG. 14. (Color online) Simulation results for the susceptibility X 2 
of model X obtained with different sample sizes; same meaning of 
symbols as in FIG. 4. 


all sample sizes to the susceptibility given by Eq. (19), we end 
up with the results of Table IV with an estimate of the transi¬ 
tion temperature &bkt = 0.275 ± 0.005 for model X. Here 
again we observe a discrepancy between the result obtained 
by fitting the bulk expression of the susceptibility to the data 
for the largest size and the ESS analysis. This may be traced 
back to the huge finite-size effects. 


TABLE IV. Estimates of 77 for Model X obtained via a log-log fit 
according to Eq. (19) for different temperatures along with the cor¬ 
responding error 5rj. 


T 

0.260 

0.265 

0.270 

0.275 

0.280 

0.285 

0.290 

V 

0.237 

0.243 

0.248 

0.250 

0.258 

0.266 

0.269 

5r] 

0.006 

0.005 

0.005 

0.006 

0.004 

0.006 

0.003 
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C. Comparisons with the regular counterparts 

As for the regular counterparts [Eqs. (4)], the existence of 
a BKT transition is by now a well-known result for planar ro¬ 
tators [Eq. (4a)], and an estimate of the transition temperature 
to be found in the Literature is Tbkt = 0.8929 ± 0.0001 
[51, 52]; Tbkt found for the V model is about 2% higher 
than the corresponding value for the regular counterpart. 

On the other hand, available evidence does not seem to sup¬ 
port a BKT scenario for the classical 0(3) Heisenberg regular 
counterpart [Eq. (4b)]. Various authors (see, e.g.. Ref. [53]) 
have argued that the model does not exhibit such a transition; 
the opposite view has been put forward by Patrascioiu and 
Seiler, in a series of papers (see e.g. [54]); examples of the 
resulting debate can be found in or via Refs. [55]. 

The nematic case [Eq. (4c)] has been studied for some 30 
years [56-67], and a BKT scenario has been proposed by var¬ 
ious Authors: a recent estimate of the transition temperature 
is Tbkt = 0.548 ± 0.002 [67], with the C* maximum at 
T" Ri 0.57, and Tct ~ T"-, on the other hand, some other 
Authors claim that the named model [Eq. (4c)] does not ex¬ 
hibit any critical transition, but its low-temperature behavior 
is rather characterized by a crossover from a disordered phase 
to an ordered phase at zero temperature [68, 69]. 


TABLE V. A summary of characteristic temperatures for the three 
models examined hy simulation in the present work; see text for def¬ 
initions. 


Model 

T' 

rjn// 

Tch 

Tbkt 

V 

« 0.75 

« 1.25 

« 1.3 

0.910 ±0.005 

W 

« 0.4 

« 0.62 

« 0.56 

0.275 ± 0.005 

X 

« 0.3 

« 0.5 

« 0.45 

0.275 ± 0.005 


These comparisons (see also Table V) suggest that, on the 
one hand, the singular character of the interaction may bring 
about a BKT behavior where the regular counterpart does not 
support it (W model); on the other hand, the effect on Tbkt 
appears to be milder where the regular counterparts already 
support this transitional behavior, and this we interpret as a re¬ 
flection of the potential features pointed out previously (Sect. 
11), in the discussion of Eqs. (5) and of EIGS. 1, 2 and 3. 

In contrast to the regular counterparts, where the temper¬ 
ature dependence of C* shows a simple maximum, upon in¬ 
creasing temperature from T = 0, the three singular models 
investigated here exhibit first a minimum and then a maximum 
of C *; this behavior also appears connected with the potential 
features discussed in Sect. II. 


V. SUMMARY AND CONCLUSIONS 

We have revisited and generalized a previously studied 
model [25, 26] and defined a few others, whose pairwise in¬ 
teractions are isotropic in spin space and restricted to near¬ 
est neighbours; in contrast to other extensively studied mod¬ 
els, their functional forms contains logarithmic singularities 


which, so to speak, do not disturb the thermodynamics. When 
D = 1, the above models could be solved in closed form, in 
terms of Gamma, Beta and Polygamma functions, and were 
found to produce orientational disorder and no phase transi¬ 
tion, at all finite temperatures, in the thermodynamic limit. 
Some of the above models have been studied by simulation 
for D = 2: among a few candidates (see Section II), we had 
chosen those functional forms which strongly favour mutual 
parallel orientations, thus strengthening (at least) short-range 
correlations; in the absence of more stringent rigorous results, 
the obtained simulation results point to orientational disorder 
at all finite temperatures, and suggest a BKT scenario in the 
three cases; we hope to carry out a more thorough simulation 
study of the models. 

Moreover, the investigated models contain logarithmic sin¬ 
gularities, causing them to slowly diverge as r ——1 or 
T —0; on the other hand, comparison with the regular coun¬ 
terparts and with the above constrained models (Section II) 
leads one to speculate as to what happens if the interaction po¬ 
tential is chosen to be more confining, i.e. made more rapidly 
divergent as t moves away from (actually, a multitude 
of such functional forms can be envisaged); preliminary work 
along these lines has been started, and its results will be re¬ 
ported in due course. 
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Appendix A: Exact solutions for D = 1 

Some available exact results in one dimension are recalled 
here; when D = 1 (hence Xj = j G Z), for a linear sample 
consisting of N spins, the Hamiltonian reads 

N 

77 = ^F(wj-Wj+i), (Al) 

i=i 

where we assume periodic boundary conditions i.e. wjv+i = 
wi; the corresponding overall partition functions can be cal¬ 
culated exactly, and this is usually realized based on the under¬ 
lying 0(n) symmetry, by means of an appropriate coordinate 
transformation (i.e., geometrically, by taking each spin Wj as 
defining the reference axis for the next one w^+i) [70-75]; the 
corresponding overall partition function reduces to the TV—th 
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power [or (iV — 1)—th power if one uses free boundary condi¬ 
tions] of a single-particle quantity, to be denoted here by g(T); 
in formulae 


,iT) = ^pin 

c2tz 


(A2a) 


p{T) = / exp(—/3F(cos 5))c?s, n = 2, (A2b) 

Jo 


and 


q{T) = ^piT), 


(A3a) 

p{T) = J exp{—f3F{s))ds, n = 3, (A3b) 

where /3 = 1/T; correlation functions are defined by 

Gj{rn) = {£j{wj ■ Wfe)), as function of m = |xj — x^j; 

(A4) 

here J is a strictly positive integer, and £j{t) denote appro¬ 
priate orthogonal polynomials, i.e. 


£j{t) = 


Tj{t) = cos(J arccos(T)), n = 2 
Pj{t), n = 3 


(A5) 


here Tj{...) denote Chebyshev polynomials of the first kind, 
and Pj{...) denote Legendre polynomials. For general D, 
and when F{t) is not an even function of its argument, the 
simplest correlation function is G'i(r); for D = 1, the defini¬ 
tion in Eq. (A4) simplify to 

Gj(m) = (f j(wj-Wfe)), as function of m = |j —fc|; (A6) 
and Gi (m) reduces to the m—th power of the quantity 

ri{T) 


ci{T) = 


P{T) ’ 


(A7) 


where 


ri(T) = / cos s exp(—/3F"(cos s))c?s, n = 2, (A8a) 
Jo 

ri{T) = J sexp(—/3F(s))c?, 


n = 3. 


(A8b) 


The corresponding susceptibility is given by [38, 39] [see also 
the following Eqs. (13c) and (14a)] 


Xi = 


I 

N 


N N 


• Wfc) 


\i=i 
N N 




N 

j=l k=l 
o N N 

j=l k=l 

hence, in the large-limit, 

, 1 + Cl 


Xi= fJ 


1 - Cl ' 


(A9) 


(AlO) 


These quantities have been calculated in the Literature in a 
few cases, where F{t) is a simple polynomial of its argument, 
i.e. F = ±T {n = 2,3), F = ±P 2 {t) (n = 3) [70-76]; in 
the latter cases F(t) is an even function of its argument, so 
that the simplest relevant correlation function is 

G 2 (m) = {P 2 (wj ■ Wfe)), as function of m = |j — fc|, 

(All) 

which similarly reduces to the m—th power of 
/ , r’2(T) 

r 2 (T) = J P2{s)exp{-/3F{s))ds, n = 3; 

(A 12b) 

in the large-limit, the corresponding susceptibility reads 
X2 = P-. -■ (A13) 

1 - C2 

Notice that the continuity of F{t) implies convergence and 
regularity of g(T); moreover the definitions entail |ci(T)| < 1 
or |c 2 (T)| < 1 at all finite temperatures; thus leading to the 
well known results related to the absence of phase transitions 
at all finite temperatures, orientational disorder in the ther¬ 
modynamic limit at all finite temperatures, and exponential 
decay with distance for the absolute value of the correlation 
functions; actually, these results may also hold under weaker 
conditions on Fir). 

There also exist in the literature a few lattice-spin mod¬ 
els involving mild integrable singularities, i.e. defined by 
bounded and generally continuous functions of the scalar 
products, which still allow usage of the method outlined here 
when 79 = 1; one such case is the sign or step model [77-83], 
defined by 

F{t) = ± sign(r); (A14) 


the model was solved exactly for D = 1 and n > 2 [79], and 
proven to remain orientationally disordered even at T = 0, 
where calculations in Ref. [79] yield for the ferromagnetic 
case 


1 r(n/2) 

" 7 ^r((n + l)/2)’ 


(A15) 


for D = n = 2 there is consistent evidence of orientational 
disorder at all temperatures, as well as of the existence of a 
BKT transition [81-83]. 

We notice in passing that other extensions of Eq. (A14) can 
be envisaged, e. g. 


F{t) = ±sign(Pj(r)), n = 3, (A16) 


where, say, J = 2,3,4; when D = 1, the resulting partition 
functions can be worked out in closed form as well. 

The effect of divergences in F{t) was seldom investigated, 
and we shall be considering here some extensions of Eqs. (3a) 
and (3b), in addition to Eq. (3c), 

Vi{t) =-ln{l +It), n = 2, (A17a) 

Wi{t) = - ln(l -f It), n = 3, (A17b) 
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where / = ±1 defines the ferro- or antiferro-magnetic charac¬ 
ter of the interaction. Both V/(r) and Wi{t) attain their min¬ 
imum when T = I, and slowly diverge to -|-oo as t — —I\ 
X{t) attains its minima when r = ±1 and slowly diverges to 
-foo as r —0; the above functions are bounded from below, 
continuous almost everywhere, and possess integrable singu¬ 
larities; moreover, their functional forms turn out to be com¬ 
putationally convenient for D = \. Two other related models 
can be defined as well, by combinining ferro- and antiferro¬ 
magnetic cases of Vj (r) with equal positive weights, and sim¬ 
ilarly for W7(r); in formulae: 

^ 2 ('r) = - ln(2(l - r^)), n = 2, (A18a) 

A 3 (t) = — ln(l — r^), n = 3. (A18b) 

Both A 2 {t) and A^(t) are even functions of their argument, 
attaining their minimum for r = 0 and diverging to -|-oo for 
|r| —-f 1; the letter A in the names recalls their antinematic 
character. Actually, further generalizations of the Vi models 
are possible, i.e. 

Vi,k{t) =+Icos{KXjk)), n = 2 (A19) 

where K is an arbitrary, strictly positive, integer, and Vf i = 
V]. By now it has been known for some time that interaction 
models only differing in the value of K produce the same par¬ 
tition functions, and that the resulting orientational properties 
can be defined in a way independent of K [21, 84, 85]; for 
more details see Appendix B. A few specific cases are listed 
here 


V+ia{t) = -ln(l -f cos(Aj;c)), (A20a) 

1Ai,i(t) = - ln(l - cos(Ajfc)), (A20b) 

V"+i, 2 (t) = - ln(l -I- cos(2Ajfc)), (A20c) 

V-i, 2 {t) = - ln(l - cos(2Ajfc)). (A20d) 

The standard trigonometric identity 

cos(2a;) = 2 cos^ x — 1 


entails that 


14 i, 2 (r) =-ln(2cos^ Ajfc), (A20e) 

^-i.2(t) = - ln(2sin2 A,fe); (A20f) 

one recognizes that V+ 1,2 defines the 2—component counter¬ 
part of the X model, and that V- 1^2 essentially coincides with 

A 2 - 

The above models can be solved explicitly, as worked out 
in the following: notice also that some qualitative results can 
be obtained in a more direct and elementary way, e.g., for 

yi.i(T-)’ 


p27r 

p(T) = / (1 +/C0S5)^(i5 

^0 

^+7r/2 

= / [(1-f coss)^-f (1 — coss)^] ds; 

J- 1^12 


(A21a) 


and, for the correlation function, 

^271- 


ri(T) = / COS s{l1 COS s)^ds 

Jo 

r-+T^/2 

/ cos s [(1 + cos sY — (1 — cos s)^] ds; 

d-7r/2 


= I 


(A21b) 

since | coss| < 1, the above equations entail |ci(T)| < 1 at 
all finite temperatures. A similar approach can be used Wj (r), 
i.e. 

p{T)= / (l-f/s)^ds 
^+1 

= / [(1 -I- sY + (1 - s)^] ds; (A22a) 

Jo 

and, for the correlation function, 

r+i 

ri{T) = / s(l-|-/s)^ds 

r + l 

= I s [(1 -f - (1 - sY] ds; (A22b) 
Jo 

since |s| < 1, the above equations entail |ci(T)| < 1 at all 
finite temperatures. 

Notice that, for each of the two functional forms (A17a) or 
(A17b), and in the absence of an external field, the two pos¬ 
sible choices for I define models producing the same partition 
functions and correlation functions related by appropriate nu¬ 
merical factors (equivalent by spin-flip symmetry). 

The above models can be solved explicitly in terms of 
known special functions with well defined analytic proper¬ 
ties, and some of them yield results involving the functions: 
Gamma 

^+00 


r 

r(z) = / 

Jo 


s^ ^ exp(—s)ds, 


Beta 


B{x,y) = B{y,x) = 


nYny) . 

r{x + y) ’ 


and Polygamma 




d'+i 


■Inr(z). 


dz'-+'^ 

Here x, y, z are complex variables with 5R(x) > 0, > 

0, 5P(2:) > 0, and I denotes a nonnegative integer [86, 87]; let 
us also recall that T (i) = y/ir. 

The above properties of V models read 

^271- 

/o 

p 2 -K 

= 2^ / (sin^ s)^ds 
^0 


p27T 

p{T) = 2^ / (cos^ sYda 
Jo 


= 2v^2'5 


np+\) 

r(/3 + i) 

q(T) = 

’ TT r(/3 + i)’ 


(A23a) 


(A23b) 
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the configurational specific heat (in units ks per particle) can 
be obtained via the appropriate derivatives of the partition 
function and reads 


C* 


1 

2^2 




(A23c) 


connected by appropriate sign factors; thus the four named in¬ 
teraction models [Eqs. (A20)] produce one and the same par¬ 
tition function, and essentially the same orientational proper¬ 
ties. 

The corresponding properties for A 3 model can be obtained 
in closed form as well; 


For V+i ,2 


^ ^ + 1 ’ 

(A23d) 

and in general for Vj^k 



(A23e) 

notice that C 2 for V+i _2 is the same as Ci for V+ 14 . 
The corresponding results for Wj (r) are 


20 

(A24a) 

“ 1 + T 2 ’ 

(A24b) 


(A24c) 

For X (r) one finds 



(A25a) 

- 1 + ^ 2 ’ 

(A25b) 


X{t) is an even function of its argument, and the previous 
Eqs. (A12a) and (A12b) specialize to 

«(!■) = (A25c) 

^ + 1 

r 2 (r) = J P 2 {s)\sfds, (A25d) 

and eventually 

(^25e) 

c,(T) = (A25f) 

Notice also that both Wi{t) and X{t) yield the same expres¬ 
sion for the configurational contribution to the specific heat 
per particle [Eqs. (A24b) and (A25b)], and produce rather 
similar expressions for ci [Eq. (A24c)] and C 2 [Eq. (A25f)], 
respectively. As for the four V models in Eqs. (A20), let us 
recall that models with the same I and different K produce 
the same partition functions, and their orientational properties 
can be defined in a way independent of K, i.e. Gmir) for 
f /,1 is the same as G 2 m{'r) for V /_2 [21, 84, 85]; on the other 
hand, the above calculations also show that V+ 1,1 and V-i^i 
produce the same partition functions and correlation functions 


q{T) 


Vtt r(^ -f 1) 

2 ^ r(/3 + |) 


(A26a) 


G* 


1 


4/1,1 + 


T 


-4-1 


’ 2 



(A26b) 


‘“O'' 

Notice that one can combine the potential models X and A 3 
to define 


r(T) =-ln[T^(l - r^)], n = 3; (A27) 


in this case the interaction diverges to +00 when r = 0 and 
|r| = 1 ; on the other hand, by standard trigonometric identi¬ 
ties, one can recognize that the n = 2 counterpart corresponds 
to 14-14 within numerical factors. The partition function of 
model Y is 


qiT) 


v/^^_g r(2/3+l) 

2 r(2/3+|)’ 


and the corresponding quantities are given by 


(A28a) 


G* 


4 

J '2 


4- 





(A28b) 


C 2 (T) 


1 

4 + 3T’ 


(A28c) 


In all of the above cases, G* was found to be a monotonic 
decreasing function of temperature, in contrast to the regular 
counterparts Eqs. (4a) and (4c), which produce a maximum of 
G* (r); on the other hand, Eq. (4b) also produces a monotonic 
decreasing behavior for C*{T). 

In the main text we are simply referring to V+ 1^1 as V 
model, and to IT+i as W model. For D = 1, the named 
models produce no phase transition and no orientational or¬ 
der at finite temperatures in the thermodynamic limit; actu¬ 
ally, some non-integrable singularities in F{t) can produce 
the same qualitative behavior as well; this happens, for exam¬ 
ple, with constrained models, defined as follows: let sq denote 
a real number, 0 < sq < tt, tq = cos sq, and let [9, 88-90] 


f fir) , + 1 > r > To 
F(t) = I , (A29) 

[ +00 , — 1 < r < To 

where / (r) denotes some regular function of its argument (see 
also below); in other words, the absolute value of the angle be¬ 
tween the two interacting unit vectors, defined modulo 27r, is 
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constrained to remain below the threshold sq. Upon follow¬ 
ing the previous line of thought and applying Eqs. (A2) to 
(A8b), one can recognize that, when D = 1, functional forms 
like Eq. (A29) also produce no phase transition and no ori¬ 
entational order at finite temperatures in the thermodynamic 
limit. Models defined by Eq. (A29) and D = n = 2 have also 
been addressed: for fir) = —r, it was proven that, when sq 
is sufficiently small, the correlation function G'i(r) never de¬ 
cays exponentially with distance, but obeys an inverse-square 
lower bound at all temperatures [9, 88, 89]; on the other hand, 
when /(r) = 0 [90], the system is athermal, and there is a 
simulation evidence of a BKT transition with sq as control 
parameter. 

Appendix B: Mapping between potential models 


where K is an arbitrary non-zero integer, and rki = -0; one 
can immediately verify that 

Vi^eZ\{O},^'K = 0; (B3) 

consider now 

p2-K 

5= / exp(±i^s)$(cos ATs, sin Ars)(is, (B4) 

Jo 

where n > 1 denotes an arbitrary positive integer, and recall 
the identity 


Consider the integral 


0 



<i)(coss, sin s)(is, 


(Bl) 


where $ denotes a sufficiently regular function, and let 


r2Tr 


'i/K = 


$(cos Ks, smKs)ds, 


(B2) 


E 

f=i 


exp ^±27ri— ) = 0. ^ > 1. 


(B5) 


Thus the value of S in Eq. (B4) is zero when p is not an in¬ 
teger multiple of K-, on the other hand, when p is an integer 
multiple of K, say /i = XK, the value of S is again indepen¬ 
dent of K 
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